Draft version June 15, 2011 

Preprint typeset using WT^^ style emulateapj v. 11/10/09 



TWO UPPER LIMITS ON THE ROSSITER-MCLAUGHLIN EFFECT, WITH DIFFERING IMPLICATIONS: 
WASP-1 HAS A HIGH OBLIQUITY AND WASP-2 IS INDETERMINATE * 

Simon Albrecht\ Joshua N. Winn\ John Asher Johnson^, R. Paul Butler^, Jeffrey D. Crane'', Stephen A. 
Shectman", Ian B. Thompson^, Norio Narita^, Bun'ei Sato'', Teruyuki Hirano Keigo Enya**, Debra Fischer^ 

Draft version June 15, 2011 

ABSTRACT 

We present precise radial- velocity measurements of WASP-1 and WASP-2 throughout transits of 
their giant planets. Our goal was to detect the Rossiter-McLaughlin (RM) effect, the anomalous radial 
velocity observed during eclipses of rotating stars, which can be used to study the obliquities of planet- 
hosting stars. For WASP-1 a weak signal of a prograde orbit was detected with w2(t confidence, and 
for WASP-2 no signal was detected. The resulting upper bounds on the RM amplitude have different 
implications for these two systems, because of the contrasting transit geometries and the stellar types. 
Because WASP-1 is an F7V star, and such stars are typically rapid rotators, the most probable reason 
for the suppression of the RM effect is that the star is viewed nearly pole-on. This implies the WASP-1 
star has a high obliquity with respect to the edge-on planetary orbit. Because WASP-2 is a KIV star, 
and is expected to be a slow rotator, no firm conclusion can be drawn about the stellar obliquity. 
Our data and our analysis contradict an earlier claim that WASP-2b has a retrograde orbit, thereby 
revoking this system's status as an exception to the pattern that cool stars have low obliquities. 
Subject headings: techniques: spectroscopic - stars: rotation - planetary systems ~ planets and satel- 
lites: formation - planet-star interactions 



1. INTRODUCTION 

The existence of Jupiter-sized planets on very close-in 
orbits presents a challenge to any model which aims to 
explain the formation of planets. In the current picture, 
these planets form further away from their host star and 
migrate inward. How and why this migratio n occurs is 
subje ct to debate (e.g. Lin et al^[T99 6; Nagasawa et al. 
2008 1 . Recently an important clue to this riddle was re- 
vealed: a subset of the close-in planets have orbits that 
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are seeming randomly-oriented with respect to the equa- 
torial plane of the host star (see, e.g. , [Hcbrard ct al. 2008 
Winn et al.|2009[| NarIta et al.|2009[ [Johnson et al.,200!^ 



[Triaud et a.. 



|2i009l|N 



' Winn et aTT l pOlOl and |Schlaufman| p)10| ) found 
that planets orbiting stars with effective temperatures 
^ 6250 K (i.e., mass ^ I.2M0) tend to have an orbital 
axis misaligned with respect to the stellar spin axis, i.e., 
a high stellar obliquity. In contrast, the two axes are gen- 
erally well-aligned for systems in which the host star is 
cooler (i.e., less massive). These authors noted that this 
could reflect a difference in the dominant planet migra- 
tion m echanism between lo w-mass stars and high-mass 
stars. Winn et al. (2010) further speculated that all 
close-in giant planets are transported inwards by pro- 
cesses that disrupt spin-orbit alignment. Subsequently, 
the angular momenta are realigned via tidal interaction, 
and this process is more rapid in cooler stars perhaps due 
to their thicker convective envelopes. In this picture any 
viable migration process would have to introduce mis- 
alignment between orbital and stellar spin. 

However, the small sample of accurate and precise 
measurements of stellar obliquities (sa25 systems) and 
the possibility of selection effects present us with many 
pitfalls if we want to validate or reject theories of giant 
planet migration. Here we report on our attempts 
to measure the spin-orbit angles in the WASP-1 and 
WASP-2 systems, taking advantage of the Rossiter- 
McLaughlin (RM) effect. 



WASP-lb was discovered by Cameron et al. (2007). 



It orbits on a 2^52 circular orbit around a F7V star and 
has a mass of 0.92 Mj„r,. One reason why this system 
is interesting is that Stempels et al. (20071 reported a 
projected stellar rotation speed of ■ysini^ < 5.79 ± 0.35 
kms^^, which is relatively slow for a star of this spectral 
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Fig. 1. — Geometry of the Rossiter-McLaughlin effect. The left panel illustrates a transit, with the planet crossing from left to 
right. Due to stellar rotation the left side of the star is moving towards the observer and the right side is receding. The unit vectors and 
fio point along the sky-projected stellar rotation axis and planetary orbital axis. They are separated by an angle A. In this diagram, 
points in the j;-direction, and the anomalous radial- velocity caused by the planet is proportional to x (see Section 2). The extrema in the 
RM signal occur at ingress {x = xi) and egress (x = X2). The relations between xi, X2, A and the impact parameter b are indicated on the 
diagram. The right panel shows the corresponding RM signal as a function of time, for an idealized case with no stellar limb darkening. 



type. For this reason, Schlaufman (2010) identified 
WASP-1 as a hkely case ot spin-orbit misaHgnment 
along the hne of sight, i.e., sini^, < 1 even though 
sinio w 1 for the planetary orbit. The star's effective 
temperature places it right in the range where the 
transition f rom well-aligned t o mi saligned orb i ts wa s 



observed Winn et al. 



Recently Simpson et al. 



piO| and |Schlaufman| ( [2010| . 

^ |2011[ ) reported a detection of 

the RM ettect tor this system and concluded the orbital 
and stellar spins were misaligned in the plane of the 
sky. As we will discuss in Section [3) our analysis leads 
to a more complex conclusion: while we agree that the 
spin and orbital vectors are misaligned, the evidence 
for a sky-plane misalignment is much weaker than the 
evidence for a line-of-sight misalignment. 

WASP-2b was also discovered by [Cameron et al.| 



(2007). This 0.87 Mjup planet has a host star of later 
spectral type (KIV) and o r bits o n a circular 2^15 or- 
bit. Recently Triaud et al. (2010) reported an angle of 
153^15 de grees between the projected orbital and stellar 
spins, i.e., a retrograde orbit. This is interesting as the 
host star is firmly on the "cool" side of the proposed di- 
vide between cool well-aligned stars and hot misaligned 
stars. WASP-2 would therefore constitute an important 
exception to the trend. However, as we will discuss in 
Section |4j we find no evidence for a retrograde orbit and 
argue that the obliquity of the host star cannot be deter- 
mined from either the new data or the previously pub- 
lished data. 

2. ROSSITER-MCLAUGHLIN EFFECT 

From the perspective of this study there are two main 
differences between the WASP-1 and WASP-2 systems. 
First, the stars are of differing spectral type, leading 
to different a priori expectations for the stellar rotation 
speed. The implications of this difference are discussed in 
Sections |3] and |4] Second, the planets' trajectories across 
the stellar disk nave different impact parameters: WASP- 
Ib nearly crosses the center of the disk, while the transit 
of WASP-2b is off-center. This section is concerned with 
the implications of this geometrical difference, as well as 



the more general relation between the characteristics of 
the RM signal and the parameters that are often used to 
model the signal. Some of these aspects of RM modeling 



were described by Gaudi & Winn 



(2007), to which we 



refer the reader for a more comprehensive account. 

Models of the RM effect wit h varying degrees of accu- 
racy h ave been worked out by IHosokawa' (1953) ; 'Quelozj 
et all (2000); Ohta et al. (2005); Winn et al. (2005) 
UmTcncz (2006); Albrecht et al. (2007); Gaudi k Wiiml 
( 12007 ); Collier Cameron et al. (2010 ) ; Hirano et al . (2010l 



and Shporer &: Brown 7 2011). Because our aim in this 
section is pedagogical, we ignore the influence of stellar 
limb-darkening, differential rotation, gravity darkening, 
surface velocity fields and any departures from sphericity 
of the planet or star. We also assume that the planet-to- 
star radius ratio Rp/R^, is small, and that this parameter 
is known precisely along with all the other parameters 
that are derived from photometric observations of tran- 
sits. In particular we assume precise knowledge of the 
impact parameter b = rt cos io/ R*, where rt is the or- 
bital distance at the time of transit, i?^ is the stellar 
radius, and io is the orbital inclination. 

With these approximations, the anomalous radial ve- 
locity due to the RM effect is 



AFrm(<) ~ - 



Rp 
R~. 



(1) 



where Vp{t) is the "subplanet" radial velocity, i.e., the ra- 
dial component of the rotational velocity of the portion 
of the photosphere hidden by the planet. Neglecting dif- 
ferential rotation, we may write 



Vp{t) = (usini^) x/R^, 



(2) 



where x is the distance on the sky plane from the center 
of the planet to t he ste llar rotation axis [see, e.g., pages 
461-462 of Gra3^(|2005|]. 

The situation is illustrated in Figure [T] In this di- 
agram, and no are unit vectors paraTTel to the sky 
projections of the stellar and orbital angular momenta, 
respectively. The angle A is measured from n^, to Uoj^ 

This definition of A is taken from |Ohta et al.| | |2005| . Some 
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Fig. 2. — The dependence of the RM signal on A, for high 
and low impact parameters. The upper left panel shows the 
geometry for a system with b 0.6, for two different cases of A. In 
the first case (solid line) the orbital and stellar spins are aligned, 
and in the second case (dashed line) they are misaligned. The 
upper right panel shows the corresponding RM signals; both the 
mean amplitude and the asymmetry of the RM signal are different. 
The two lower panels show a similar orbital configurations but 
for 6 0. Here, the mean amplitude changes with A but the 
asymmetry is always zero. 



Fig. 3.— Spectrum of the F7V star WASP-1. A small 
portion of the the spectrum of WASP-1, as obtained with HIRES, 
is shown. The dots represent the the observed spectrum, the solid 
line represents our best fit with a macro-turbulence parameter of 
3.98 kms~^ and a iisini* of 2.9 km s^^. The (red) d ashed line 
shows the s pectrum broadened to the values given by [Stempels] 
|et al. | |2007 1, who obtained a usini^, of 5.79 km s~^ with a macro- 
turbulence parameter of 4.5 kms~^. The two lower rows of points 
show the differences between model a nd data f or our best fit (black 
dots) and the values given by jStempels et aL| | |2007| (red dots). 



The maximum redshift and blueshift occur at ingress and 
egress, which we take to have x-coordinates of xi and X2 
respectively. Using the geometrical relations shown in 
the diagram, we may write xi and X2 in terms of b and 
A, 

xi = ^\/l — — 6 tan A^ cos A — \/l — cos A — 6 sin A, 
= (Vi - + &tanA) cosA = \/l - b^ cos A + bsinA. (3) 



X2 ■- 



It is instructive to examine the (scaled) sum and differ- 
ence of xi and X2, 

- V sin {x2 + xi) = \f\--\P- v sin cos A, 

- V sin ii, {x2 — xi)~b V sin i^, sin A. (4) 

The sum is the mean amplitude of the red and blue peaks 
of the RM effect, while the difference is a measure of 
asymmetry between the peaks. For a fixed 6, the mean 
amplitude depends on v sin cos A while the asymmetry 
depends on v sin ii, sin A. 

Figure [5] shows the RM signal in 4 different situations: 
two different values of A for each of two different impact 
parameters. The upper panels show the case h ~ 0.6, as 
is the case for WASP-2. Here, as A is varied, both the 
mean amplitude and asymmetry of the RM signal are 
observed to change. By measuring the mean amplitude 
and asymmetry, one may determine both vsinii, and A. 
The lower panels show the case 6 w 0, as is the case for 
WASP-1. Here, the asymmetry vanishes regardless of 
A. The only observable quantity is the mean amplitude, 
and therefore the only parameter combination that can 
be determined is v sin cos A. 

other investigators measure the angle from no to and denote 
the angle (5. Clearly /3 = —A. 



Consequently, for transits with low impact parameters, 
A and v sin i^, have strongly correlated uncertainties and 
it is not possible to measure A without some prior infor- 
mation about vsinii,. However, in such cases it is still 
possible to tell whether cos A is positive or negative, and 
therefore whether the orbit is prograde (|A| < 90°) or 
retrograde (|A| > 90°). We also note that the degen- 
eracy between tisini^, and A can be broken in principle 
when the RM effect is modeled at the level of spectral-line 
distortion, ra ther than model ing onl y the anomalous ra- 
dial velocity (Albre cht et al.|2007[ [Collier Cameron et al. 
20101. In this paper, though, we work with the anoma- 
lous radial velocity. 

3. WASP-1 

3.1. Observations and basic stellar parameters 

We conducted spectroscopic observations of WASP-1 
transits with the Keck I 10 m telescope and the Subaru 
8.2 m telescope. With Kec k, we used t he High Resolu- 
tion Spectrograph (HIRES; |Vogt et ain[r994j to gather 
34 spectra spanning the transit of 20U7 September 1/2. 
With Subaru, we used the High Dispersion Spectrograph 
(HDS; Noguchi et al.|2002 ) to observe two different tran- 
sits, on the nights of 2007 August 4/5 and 2007 Septem- 
ber 6/7. A total of 23 spectra were obtained with HDS, 
most of which (20) were obtained on the latter night. 
At both observatories an iodine gas absorption cell was 
used to correct for changes in the point spread function 
and wavelength scale. Radial velocities (RV) were de- 
rived from t he spectra using pro cedur es similar to those 
desc ribed by [B utler et al. (1996). See [Sato et aL] ( |2002 | 
and Narita et al. (2007) for details ori the Subaru data 
reduction, i'he RVs are shown in Figure [5] and given in 
Table l 

To check on the basic stellar parameters, we used 
the Spectroscopy M ade Easy (SME) software package 
( Valenti fc PiskunoV|1996i ) to model the high-resolution. 
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Fig. 4. — Photometry of WASP-1 transits. The upper 
panel is a composite z'-band ligli t curve based on our data and 
that of [Charbonneau et al.| | |2007| . The lower three panels show 
the residuals between eacn ot tne 3 datasets and the best-fitting 
model. 




time [tir] 

Fig. 5. — Spectroscopy of WASP-1 transits. The radial 
velocities measured before, during, and after transit are plotted as 
a function of time from inferior conjunction. Solid symbols are data 
from HIRES and open symbols arc data from HDS. Gray symbols 
are the SOPHIE data from Simpson ct al. (2011), which are shown 
for comparison only (they wore not used in our fitting process). The 
upper panel shows the measured RVs and the best-fitting model. In 
the middle panel, the orbital contribution to the observed RVs has 
been subtracted, isolating the RM effect. The lower panel shows 
the SOPHIE RVs after subtracting our best-fitting orbital model. 
The light and dark gray bars in the two lower panels indicate times 
of first, second, third, and fourth contact. 



high-signal-to-noise ratio template spectrum. We ob- 
tained Teff = 6213 ± 51 K, logg = 4.19 ± 0.07, [M/H] = 



0.17 ± 0.05, and vsini^ = 1.60 ± 0.50 km s^ 



These 



can be comp ared to the previous sp ectroscopic results of 
WASP-1 by [Stempels etal] ( |2007[ ), which gave Teff = 
6110 ± 45 K logs = 4::2S±\n^M/H] = 0.23 ± 0.08, 
and vsmi^^ = 5.79 ± 0.35 km s~^. Our analysis gave a 
higher value of Toff and a lower value of v sin i^, . 

The discrepancy in T^s is discussed in Section [5] The 
discrepancy in v sin j^, is of immediate importance be- 
cause stellar rotation is a key parameter in the inter- 
pretation of the RM effect. Frequently, such discrep- 
ancies arise because of differing assumptions regarding 
turbulent broadening. SME determines wsini^, based on 
the observed widths of numerous weak lines in the spec- 
trum. The widths are influenced not only by rotation, 
but also by random motions of the stellar photosphere 
(microturbulence and macroturbulence) , and these ef- 
fects cannot generally be disentangled. Hence it is nec- 
essary to assume "typical" values of the turbulence pa- 
rameters and attribute the excess broadening of the ob- 
served lines to rotation. When using SME, it is assumed 
I'mic = 0.85 km s~^ and 

one , Teff- 5770 _i 
+ 650 K ""^^ ' 



an empirical relation determined by |Valenti &: Fischer] 
( |2005[ )P] For WASP-1, this formula gives w^ac = 



^^_Thc equation given here corrects a sign error in Equation (1) 
of |Valenti fc Fischer^ (2005J . 



4.66 km s 



This is not too different from the value 



4.5 km s that was assumed by Stempels et al. 



( 2007 1 and hence the discrepancy in v sinii, cannot be 
attributed to different assumptions regarding macrotur- 
bulence. 

To investigate further, we performed a differential as- 
say for rotation, based on a comparison between the So- 
lar spectrum and a Keck/HIRES spectrum of WASP-1. 
First, we deconvolved the WASP-1 spectrum to remove 
the instrumental broadeni ng of width 2.2 km s ~^. Then, 



using the MORPH code of Johnson et al. 
plied a rot ational broadening k ernel to 



(20061, we ap- 
hTWBO Solar 
spectrum of Kurucz et al. ( 1984 ) to achieve the best fit to 
the deconvolved WASP-1 spectrum. We found that the 
best-fitting broadening kernel was 2.36 km s^^, indicat- 
ing the WASP-1 lines are slightly broader than the Solar 
lines. Figure |3] shows a small portion of the WASP-1 
spectrum and our best-fitting model based on the broad- 
ened Solar spectrum. 

The larger breadth of the WASP-1 lines could be in- 
terpreted as more rapid rotation than the Sun, but in 
fact part of the increased breath is expected to be due to 
the higher macroturbulence of WASP-1. However since 
the accuracy of Eqn.[5]is not known, we may here simply 
assume that the macroturbulence of WASP-1 is greater 
than or equal to the macroturbulence of the Sun. The 
MORPH finding implies 

[vsim, (Wl)]2 + [t;^,, (0)]2« 

[vsim^ (0)]2 + [v^,, (0)]2 + (2.36 km s'^)^, (6) 
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Fig. 6. — Results for jjsini, and A, based on our MCMC 
analysis in the WASP-1 system. The gray scale indicates 
the posterior probability density, marginalized over all other pa- 
rameters. The contours represent the 2-D 68.3%, 95%, and 99.73% 
confidence limits. The one-dimensional marginalized distributions 
are shown on the sides of the contour plot. A strong correlation 
between the projected rotation speed and the projected angle be- 
tween the stellar and orbital spins exits. Either the two axes are 
nearly perpendicular on the sky plane, or else v sin is small and 
A can have any value. 
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Fig. 7. — Results for Dsini* and A, this time including 
a prior constraint on tjsini*. The prior constraint was based 
on the spectroscopic result vsini^, < 2.9 km s~^ (see Section 3.1) 
and is illustrated by the dashed line in the right-hand side panel. 
For wsinj* < 2.9 km s~^, the prior was set equal to unity; and 
for greater values the prior was a Gaussian function with mean 
2.9 km and standard deviation 0.5 km s~^. Compared to 
Figure |6] the solutions are similar but are constrained to have 
somewhat lower vsin^^,. 



where the "Wl" quantity is for WASP-1 and the 
"0" quantities are for the SunFj Taking the disk- 
integrated rotation and macroturbulence of the Sun to 
be 1.63 km and 3.98 km s~^, and the macroturbu- 
lence for WASP-1 the same as the sun, Equation |6] gives 
vsini^ < 2.9 km s'^ for WASP-1. 

These results show that the projected rotation speed of 
WASP-1 is quite slow (<2.9 km s~-^) and is in fact nearly 
undetectable against the dominant line-broadening effect 
of macroturbulence. Figure |3] also shows that our spec- 
trum is incompatible with the more rapid rotation of 
5.79 ± 0.35 km s~ ^ found by IStempels etall ( |2007[ ). We 
do not know why Stempels et al. 1 2(Jl)7 ) found a higher 
V sin even when making equivalent assumptions regard- 
ing macroturbulence. Genuine changes in vsini^, could 
be produced by spin precession, but are not expected to 
be appreciable on such short timescales , and hence we 
procee d under the assumption that the Stempels et al. 
(2007) determination was in error. 

lb reduce the uncertainties in the photometric pa- 
rameters we gathered new photometric data with Ke- 
plercam, a CCD camera on the 1.2 m telescope of the 
Fred L. Whipple Observatory on Mount Hopkins, Ari- 
zona ( Szentgyorgyi et ar]|2005 1 . Observations were con- 
ducted in the SUSS 2;' -band on 2009 September 16/17 
and 2010 September 29/30, although bad weather inter- 
rupted the transit in both cases. The new photometric 
data were combin ed with the previous data of |Charbon-| 
neau et al. (2007), which were gathered with the same 
instrument and reduced with similar procedures. All of 
the Keplercam data are shown in Figure |4] 

13 verified with numerical experiments that in this regime of 
velocity widths and for the SNR and resolution of our spectrum, 
the widths of the various convolution kernels can be approximately 
added in quadrature as implied here. 



3.2. Analysis 

To derive constraints on A, we fitted a model simulta- 
neously to the RV data and the photometric data . The 
photometric trans it was modeled with the code of 



Man- 



del fc Agoll ( |2002[ ), and the RM e ffect was modeled with 
a simplihed version of the code of Albrecht et al. (2007). 
This model for the RM effect is similar to that given m 
Eqn.[7]but takes limb darkening into account. It does not 
take into account the nonlinear relation between AVrm 
and Vp(t) because those nonlinearities are important only 
for stars with larg er vsmi^, (see, e.g. 
Hirano et al]|2010| ) 



Winn et al. 2005 



I'h e transit imp act parameter for WASP- lb is small. 



with 



Torres et al. 



1+0.27 
-0.00- 



mm having reported b = O.OOZ _ 
Therefore, based on the reasoning of Section[2j we expect 
the data to constrain tisini^^cosA but not w sin i^, sin A. 
For this reason we chose to parameterize the RM ef- 
fect with the quantities \/v sini^, cos A and \/v sini^, sin A, 
rather than wsini^, and A. The reason for the square 
roots is to give a constant Jacobian between the fitting 
parameters and the "physical" parameters vsinii, and A. 
As a result, uniform priors in our fitting parameters cor- 
respond to the desired uniform priors in usini^ and A. 
With no square roots, and no other adjustment to the 
fitting procedure, the implicit prior would be linear in 
vsinii, and would thereby bias the results toward faster 
rotation rates. 

The other model parameters were a constant RV offset 
specific to each spectrograph; the semiamplitude of the 
star's orbital velocity (iiT*), which controls the RV slope 
that is observed on each transit night; the orbital pe- 
riod (P); a particular time of midtransit (Tc); the stellar 
radius in units of the orbital distance (i?^,/o); the co- 
sine of the orbital inclination (cosio); the planet-to-star 
radius ratio [Rp/ Ri,)] two quadratic limb-darkening coef- 
ficients Ml and U2 for describing the z'-hand photometric 
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TABLE 1 

Parameters of the WASP-1 system 



Parameter Values 



Parameters mainly derived from photometry 



Midtransit time Tc [BJDtdb-2 400 000] 


54461.8630 ± 0.0002 


Period, P [days] 


2.5199464 ± 0.0000008 


cos io 


0.000-0.034 


Fractional stellar radius, Rir/a 


n IVQ + 0.003 


Fractional planetary radius, Rp/Ri, 


0.1059 ± 0.0006 


U\+U2 


0.20 ± 0.05 


Parameters mainly derived from RVs 


Velocity offset, HDS [ms"!] 


± 1.5 


Velocity offset, HIRES [ms^^] 


-17 ± 2 


Velocity semiamplitude, Ki, [ms"-*^] 


125 ± 5 


v'l' sin sin A [kms"-*-] 


-0.6 ± 0.9 


sin cos A [kms~^] 


0.31 ± 0.25 



Indirectly derived parameters 



Orbital inclination, io [°] 88-92 
Full duration, T14 [hr] 3.684 ± 0.017 

Ingress or egress duration, T12 [min] 21.5 it q I 

Projected stellar rotation speed, usini* [kms^^] 0.7 it q | 

Projected spin-orbit angle, A [°] —59 it 



data; and a linear limb-darkening coefficient u to describe 
the spectroscopic transit (for which most of the signal 
is derived fro m the r egion 5000-6200 A). According to 
the tables of iClaret ( 2004[ ), appropriate choices for the 
limb-darkening coetiicients are Mi=0.1666, M2=0.3583, 
u = 0.6. We allowed ui + U2 to be a free parameter 
and held fixed ui —U2 at the tabulated value of —0.1917, 
since the difference is only weakly constrained by the 
data (and in turn has little effect on the other parame- 
ters). Likewise we held fixed u ~ 0.6. We assumed the 
orbit to b e circular, as no sign o f any e ccentric ity was de 



tected by Cameron et al 



(2009), Wheatley et al. (20101, or Pont et al. (2011) 



any ( 

(120071 , [Madhusudhan & Winn 



All 01 the time stamps of the spectroscopic and photo 
metric data wer e placed on the BJDt db system using 



the algorithm of Eastman et al. 
The fitting statistic was 



(20101. 



+ 



RV,(o) -RV^c) 

- 115ms~ ^^ ^ 
llms-i 



t2 1134 

+ E 

j=i 



F,(o)-F,(c) 



(7) 



where the first two terms are sums-of-squares over the 
residuals between the observed (o) and calculated (c) val- 
ues of the radial velocity (RV) and relative flux (F) , and 
the last term re presents a prior constraint on Ki, based 
on the results of |Cameron et al.] ([2007 ) . Below we will re- 
peat the analysis including the constrain on v sin i^, found 



In particular, [Madhusudhan &: Winn| l |2009[) reported an up- 
per limit of e < 0.U8S with 95.4% confidence. If the orbit were 
actually eccentric, in contradiction of our modeling assumption, 
then the main change would be that our result for the velocity 
semiamplitude would be biased. The results for the spin-orbit 
parameters would not be significantly affected. 



in Section [XT] 

We solved for the model parameters and their uncer- 
tainties us ing the Markov Chai n Monte Carlo (MCMC) 
algorithm ( Tegmark et al.|2(504 ) . We used a chain length 
of 2 X 10^ steps and set the size of the steps in each pa- 
rameter yielding an acceptance rate of about 30%. Before 
running the chain we increased the uncertainties of the 
HIRES RVs by adding a "steflar jitter" term of 5 ms"^ 
in quadrature to the internally-estimated uncertainties. 
This choice of jitter term produced a reduced of unity 
when that data set was fitted alone. In making this step 
we have assumed that the extra RV noise is well described 
as Gaussian and uncorrelated. This is consistent with the 
appearance of the residuals shown in Figure [5j although 
we acknowledge there is no guarantee. Table [3] reports 
the original, internally-estimated uncertainties without 
any jitter term. 

The results for the RM parameters are displayed in 
Figure l6j and the results for all the parameters are given 
in Table fl] As anticipated, the weak detection (or non- 
detection) of the RM effect led to tighter bounds on 
u sin i.^ cos A than on w sin i^, sin A. This is why the con- 
tours in Figure [6] reach to large values of v sin for small 
values of cos A (A w ±90°). 

In an attempt to break the degeneracy between v sin ii, 
and A we refitted the data with a prio r co nstraint on 
wsini,^. Based on the results of Section [XTj we used a 
one-sided Gaussian prior, taking the value of unity for 
vsmii, < 2.9 km s~^ and falling off as a Gaussian func- 
tion with cr — 0.5 kms~^ for higher values. The results 
from this more constrained MCMC analysis are shown 
in Figure |6| The modified bounds on A are — 53±2g°. 
This analysis disfavors A w ±90° as this would require 
larger vsini^. However it is not possible to tell defini- 
tively whether the positive or negative solution is correct. 
Within the 95% confidence contour, all prograde orbits 
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are allowed. 

A different approach is to use a prior constraint on v, 
the actual rotatio n speed of the star , based on its spec- 
tral type and age. Schlaufman (20101 recently presented 
a formula for a mam-sequence star's expected rotation 
period, given its mass and age. He based the formula on 
the observed r otation periods of stars in young clusters 
along with the Skumanich ( 1972 ) law v oc He fur- 

ther showed that this formula gives a good description 
of the V sin ^^, d istribution of stars in the SPOCS catalog 
( [Valenti fc Fischer, 2005) . For WASP-1, he found an ex- 
pected value V = 8.6 ± 0.5 kms~^ where the uncertainty 
is based only on the uncertainties in the age and mass 
of WASP-1, and does not account for any uncertainty 
due to intrinsic scatter in the mass-age-period relation, 
which seems to be about 3 times larger than the for- 



mal uncertainty [see, e.g.. Fig. 3 of Schlaufman (20101]. 
Taking v = 8.6 ± 1.5 kms~^ together with our result 
usini^ < 2.9 kms^^, the implication is sini^, < 0.34, 
i.e., the star is viewed close to pole-on. 

One might wonder if the Skumanich law is really ap- 
plicable to stars with close-in planets, which may have 
undergone significant evolution due to tidal interactions. 
For the case of WASP-1, at least, there is supporting evi- 
dence for relatively rapid rotation, based on its observed 
color and chromospheric emission. Aigrain et al.j (|2004|) 
explain how to use a star's observed i/^T/ and logioR'^^ 
indices to predict its rotation period. Applied to WASP- 
1, for which B-V = 0.53 and logwR'n^ = -5.114 (Knut- 
son et al. 2010), we find a rotation period of 12.9 days. 
Using a stellar radius of 1.45 Rq (Charbonneau et al. 



20071, the predicted rotation speed is u = 5.7 kms ^, in 



good agreement with the value expe cted from the statis- 
tical analysis by Schlaufman (20101. 

We therefore have two inde]Dendeiit lines of evidence for 
a high obliquity, or equivalently, we have strong evidence 
against the well-aligned scenario in which sin i^, « 1 and 
A « 0°. (1) The absence of a strong RM effect re- 
quires either that |A| ~ 90°, or else usini^ is very low 
(<1 km s~^). The latter possibility is incompatible with 
a well-aligned star (sinj* ~ 1), because the rotation rate 
for a star of the given mass and age is expected to be 
8.6 ± 1.5 kms~^. The observed color and chromospheric 
activity level also suggest a rotation speed of this order. 
(2) Independently of the RM effect, our determination 
of usini^ based on the observed width of the spectral 
lines is much lower than the value of the expected ro- 
tation speed, which implies a low sinij,. In short, it is 
likely that the stellar and orbital spins are misaligned 
along the line of sight, and it is possible that they are 
also misaligned within the sky plane. 

3.3. Comparison with previous results 



Simpson et al. (2011 1 reported A=— 79 



1-1-4.5 o 
-4.3 



^ ^ ^ for WASP- 

Ib, based on observations taken during and after a plane- 
tary transit with the SOPHIE spectrograph on the 1.93m 
telescope of the Observatoire de Haute-Provence. Their 
value for A is compatible with our result. However their 
uncertainty is much smaller than we have found. What 
causes this difference in obtained confidence intervals? 

Their RV data, reproduced in the bottom panel of Fig- 
ure [5] appears to have a higher amplitude than was seen 
in our data. This could lead to a somewhat higher result 
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Fig. 8. — Spectroscopy of WASP-2 transits. Similar 
to Figure [5] Black symbols are PFS data, and open symbols are 
HDS data. Gray symbols are the HARPS data of Triau d et al.| 
| |2010| l, which are shown for comparison but were not used during 
the htting process. The upper panel shows the data and the best- 
fitting orbital model. In the lower two panels, our best-fitting 
orbital model has been subtracted from the data. 

for V sin i^ but would not by itself affect the very strong 
correlation between v sini^, and A. Rather, the important 
differences are in the methods of analysis. There are two 
main differences. 

Firstly, rather than jointly fitting the photometric and 
spectroscopic data as we have done, Simpson et al.| ( 20Tl] ) 
fitted their spectroscopic data using independent Gaus- 
sian priors on the photometric parameters a/i?^, Rp/R^,, 
and io- The problem is that those parameters are them- 
selves very strongly correlated and their posterior distri- 
butions are far from Gaussian. In particular their photo- 
metric priors excluded very low impact parameters, while 
we find that 6 « is allowed. To avoid this problem it 
is better to analyze photometric and spectroscopic data 
together, or to place priors on the relatively uncorrel ated 
parameters T14, T12 and Ru/R-i, ( Carter et al.||2008 ) . 

Secondly, Simpson et al. (20111 used a prior on vsini^, 
based on the spectroscopic analysis of |Stempels et alT] 
(|2007|, which gave wsini* = 5.79 ± 0.35 kms-\ As 



explained in Section 3.1 and shown in Figure [3] our 
spectroscopic analysis implies a slower projected rotation 
rate. Their prior on v sin pushed their solution towards 
higher v sin and excluded aligned configurations of the 
projected axes. 

4. WASP-2 

4.1. Observations and basic stellar parameters 

We conducted spectroscopic observations of WASP-2 
transits with the Magellan (Clay) 6.5 m telescope and 
the Subaru 8.2 m telescope. With M agellan we used the 



Planet Finding Spectrograph (PFS; Crane et al. 



20101 



to gather 35 spectra spanning the transit of 2010 Au- 
gust 26/27. With Subaru we used the HDS to obtain 21 
spectra spanning the transit of 2007 September 4/5, and 
10 spectra spanning the transit of 2007 September 19/20. 
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Fig. 9. — Results for I'sini* and A in the WASP-2 system. Similar to Figure|6] but for WASP-2. The gray scale plots indicate the 
posterior probability densites, marginalized over all other parameters. The contours represent the 2-D 68.3%, 95%, and 99.73% confidence 
limits. The one-dimensional marginalized distributions for A are on top of the contour plot. The left panels show the results for the MCMC 
analysis with no prior applied to Ki,, the middle panel shows the results with a prior on ftT* as shown in equation [s] and the right panel 
shows the result for the prior wit h half the confidence interval. The black dots with error bars mark the results by Triaud ct al. (2010jl. 
The error bars are those quoted by |Triaud et al.^ ( ^2010j ), representing 68.3% confidence intervals in A and tisini ^ marginafized over all other 
parameters. They are not strictly appropriate tor this two-dimensional plot. We refer the reader to Figure 3 of [Triaud et al.| | |2010| to view 
their two-dimensional posterior distribution. 



Again we employed the iodine-cell technique to derive 
precise radial velocities. All the RVs are given in Table|4] 
and plotted in Figure |8] As was the case for WASP-1, 
we found no clear evidence for the RM effect. 

To check on the basic stellar parameters, we 
also obtained a high-quality template spectrum with 
Keck/HIRES, so that we could use the same SME-based 
analysis that was used for WASP-1. Wc obtained Tcff = 
5206 ± 50 K, logf? = 4.51 ± 0.10, [M/H] = 0.04 ± 0.05, 
and usini^ = 1.3 ± 0.5 km s^^. The assumed macro- 
turbulent velocity was 3.11 km s^^. Using the MORPH 
code described in Section 3.1, we found that the WASP-2 
lines are no broader than the Solar lines, and estimate 
wsini^ ^ 1.5 kms~^. 

4.2. Analysis 

T he transit impact p arameter for WASP-2b is large, 
withlTorres et al.''''^""'''* ^--^ ^ - mA+a.on 



( 2008 ) having reported b = 0.7241o (,28- 
Theretore, based on the reasoning of Section 2, the non- 
detection of the RM effect implies that both v sin z^, cos A 
and 11 sin z^, sin A are small, which is only possible for low 
usini^. Unlike the case for WASP-lb, the RM effect for 
WASP-2 cannot be suppressed by having the planet's 
trajectory coincide with the sky-projected rotation axis. 
We therefore expect the nondetection to lead to an upper 
limit on vsinii, and no information about A. 

For the quantitative analysis our procedure was similar 
to that used for WASP-1. The RVs were modeled as 
sum of contributions from a circular orbit, the RM effect, 
and a constant offset spe cific to each spectro graph. We 
used a prior on K^, from Triaud et al. (20101, but with 
a doubled uncertainty (see below), and also tested the 
sensitivity of the results to this prior as described below. 
Since the photometric parameters are already precisely 
determined and we do not have any new photometric 
data, we implemented priors on the full transit duration 
(T14), the ingress or egress duration (712) , the radius 



ratio (Rp/R^) from Charbonneau et al. (2007 ), and the 
transit ephemeris based on the analysis of [ijouthworth 



et al. (2010). The fitting statistic was 



RV,(o) -RV.(c) 

0"RV,i 



^ ^ Tc,B.iD - 2453991.51530 y ^ /P - 2'? 15222144 A' 



0.00017 



0^00000039 / 



_^ ' T14, - 1.799hr \ ^ fTi2 - 24.6min 



0.0035 hr 



2.4 min 



+ 



Rp/R^ - 0.1309 Y ^ f K^,- 153.6 ms~^ 



0.0015 



6 ms 



(8) 



where the symbols have the same meaning as in section 
[3} For the PFS data, a "stellar jitter" term of 10 ms~^ 
was added in quadrature to the internally-estimated un- 
certainties to give a reduced of unity. This probably 
reflects the limitations of the current algorithm that is 
used to estimate uncertainties, which is geared toward 
much brighter stars. 

Our results are presented in Table [2] and are illustrated 
by the contours in the middle panel of Figure lo] [The 
single solid point in Figure [9] represents the resulx of jTri- 



aud et al.| ( |2010[ ) , which will be discussed below.] As 
expected, v sm is constrained to low values but A can 
assume any value from —180° to -1-180°. 

The three different panels of Figure |9] show the results 
of different choices for the prior on i^^*. We wondered 
about the sensitivity of the results to this prior because 
the star is a late-type star and might be expected to 
have starspots, which can cause the observed RV slope 
surrounding the transit phase to be steeper than one 
would expect from the spectroscopic orbital parameters. 
Starspots always move across the stellar disk from the 
approaching limb to the receding limb, and thereby pro- 
duce an RM-like effect with a negative slope, which is 
added to the actual orbital velocity gradient. This effect 
can be seen in a number of RM datasets presented in 
the literature, most notably for the highly spotted star 
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TABLE 2 

Parameters of the WASP-2 system 



Parameter 


Values 


Parameters mainly controlled by prior knowledge 


Midtransit time Tc [BJDtdb-2 400 000] 


53991.51530 ± 0.00017 


Period, P [days] 


2.15222144 ± 0.00000040 


cos io 


0.091 ± 0.007 


Fractional stellar radius, R^/a 


0.125 ± 0.005 


Fractional planetary radius, Rp/Rt, 


0.1309 ± 0.0015 


Parameters mainly derived from RVs 


Velocity offset, PES [ms^^] 


-2 ± 2 


Velocity offset, HIRES [ms^^] 


-23.6 ± 2 


Velocity semiamplitude, Ki, [ms^-*^] 


164 ± 4 


V'^ sin sin A [kms"-*-] 


-0.02 ± 0.28 


•^vsinii, cos A [kms"-*-] 


-0.038 ± 0.36 


Indirectly derived paramet 


ers 


Orbital inclination, io [°] 


84.8 ± 0.5 


Full duration, T14 [hr] 


1.799 ± 0.037 


Ingress or egress duration, T12 [min] 


24.2 ± 2.4 


Projected stellar rotation speed, tjsinu [kms~^] 


<0.5 (2(t) 


Projected spin-orbit angle, A [°] 


all values allowed 



CoRoT-2 d Bouchy et al]|2008[ ). Depending on the distri- 
bution of measurements betore, during and after transit 
this might introduce different biases in the resuhs for A 
and V sin i^, . 

In Figure |9j the left panel shows the results with no 
prior on if^,, the middle panel shows the result for a 
prior on Ki, as in Eqn. ^ , and the right panel employed 
the same prior but with a width of 3 m s""'^ instead of 
6 m s~^. Evidently the results are not very sensitive to 
the prior on K^,: in all cases vsinii, must be low and A 
may have any value. For concreteness our final results 
given in Table [2] are based on a prior with a width of 
6 m s~^ (i.e., the analysis depicted in the middle panel). 

One interesting feature of Figure |9] is that the poste- 
rior probability density for A has peaKS near 0° and 180°. 
For these choices of A, larger values of wsini* are com- 
patible with the nondetection. This is a general result 
when fitting RJVI data with a low signal-to-noise ratio of 
a high-& system, and can be understood as follows. For A 
near 0° and 180°, the RJVI signal is antisymmetric about 
the midtransit time. In such cases iisini^ and K^, are 
strongly correlated parameters, since small changes in 
either parameter produce changes to the RM signal that 
are antisymmetric about the midtransit time. This leads 
to larger confidence intervals for usini^. In contrast, for 
A = ±90° the KM signal is symmetric about the mid- 
transit time; it is a pure redshift or blueshift. Here, the 
parameters K^, and v sin i^, are uncorrelated and the al- 
lowed region for z;sini^, shrinks. To put it another way: 
by fitting for the systemic velocity and K^,, we have ef- 
fectively applied a high-pass filter to the RV data, and 
thereby reduced the amplitude of any RIVI signal with 
A = ±90° in comparison to the higher-frequency signal 
that is produced with A near 0° and 180°. This causes 
the allowed range of v sin to be higher for A near 0° 
and 180°. This explanation was confirmed with further 



numerical experiments described in Section 4.3 

As with WASP-1, one may try to gain more infor 
mation on the spin orbit alignment by using prior con- 
straints on vsini-i, or v, but in this case not much re- 
finement is possible. The analysis of the WASP-2 tem- 
plate spectrum gives an upper limit wsini^ ^1.5 km s^^ 
which is not constraining in this context. Also, there 
have been no reports of photometric variations due to 
star spots, and hence no stellar rotation p eriod has been 
determined. Likewise, Schlaufman (20101 found that the 
expected rotation speed tor this systein, based on its 
mass and age, is 1.61 kms~^ with an uncertainty range of 
1.72 kms~^ (presumably an asymmetric error interval). 
Because of the large uncertainty it is not possibl e to draw 
any c onclusion about sini^,, and for this reason [Schlauf-| 
man ]([2010) did not identify WASP-2 as a probable case 
ot a misaligned star. 

As an additional check on the expec ted stellar ro- 
tation speed, we used the approach of Aigrain et al 



(2004) to estimate the rotation period of WASP-2, as 
T^re^d for WASP-1. In this case, B -V = 0.84 and 
logioR'sK = —5.054 (Knutson et al. 2010), from which we 
derive a stellar rotation perio d of 46 days. Together with 
an stellar radius of 0.81 R© ( Charbonneau et al.||2007" ), 
this gives a rotation speed of v = 0.9 kms~^ which is in 



line with the low speed predicted by Schlaufman (2010). 



4.3. Comparison with previous results 

A tr ansit of WASP-2 was observed by [Triaud et al. 
( |2010[ ) with the HARPS spectrograph. Their data are 
shown in the bottom panel of Figure [8j Based on the 

-153^^^^ degrees (a retro- 



O.QQtoil kms-^ Our data are 



HARPS data they found A 
grade orbit) and v sin i^, 

not compatible with those parameters. When we fixed 
A and v sin at the values found by these researchers, 
and refitted our data, the minimum rose from 60.9 
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to 72.6, giving Ax^ = 11.7. What can have caused the 
difference between our results and theirs? 



Daemgen et al. ( 2009 1 found that WASP-2 has a neigh- 
boring star (a companion or chance alignment) at an 
angular separation of 0.7 arcsec, close enough to have 
been possibly included within the spectrograph slit or 
fiber in some cases. It is hard to predict the exact effect 
that the additional starlight would have on the spectro- 
scopic analysis, but as the neighbor is 4 mag fainter than 
WASP-2, and as its spectral type and systemic velocity 
are likely quite different from that of WASP-2, we con- 
sider it unlikely that variable contamination by this star 
is responsible for the differing results. We are therefore 
led to look elsewhere for an explanation. 

One relevant differe n ce in the analysis procedures 
is that 



Triaud et al. 



(2010) used uniform priors in 
V sin sin A and v sin i^. cos A, thereby adopting a prior 
that is linear in t>sini,^. This is in contrast to our prior 
which was uniform in vsmi^^. Their prior pushes vsinii, 
to higher values and therefore pushes A near 0° or 180° 
(see Figure [9]) . When we refitted their data using our 
procedure, we found a lower v sini^, and an enlarged con- 
fidence interval, as expected. The open circle and the 
thick dashed lines in Figure [TO] represent our fit to the 
HARPS data. However this oiiference in priors cannot 
explain the entire discrepancy: even our reanalysis of the 
HARPS data gives A = —151^^3 degrees and wsini^ = 
0.84 ±0.35 kms"^ 

This apparently statistically significant result is sur- 
prising since the RM effect is not apparent by visual in- 
spection of the data (Figure Jsl). The data during the 
transit does not appear too dmerent from the data out- 
side of the transit. If the RM effect had been measured 
but not modeled, then one would expect the residuals be- 
tween the data and the best-fitting orbital model would 
have a higher scatter inside the transit than outside the 
transit. For our data this is not the case. For the HARPS 
data set the rms residual of the out-of-transit data is 
6.9 ms~^, as compared to 7.2 ms~^ during transit. This 
represents only a marginal increase in scatter. 

This led us to conduct some numerical experiments on 
fitting random noise with similar characteristics to the 
HARPS data. We used the timestamps of the HARPS 
transit-night data, and simulated RV data based on only 
the best-fitting orbital model for WASP-2. We added 
Gaussian "measurement" uncertainties with a standard 
deviation of 7.0 ms~^. Then we fitted this mock dataset 
together with the photometric priors using a Levenberg- 
Marquardt least-squares minimization routine. This was 
repeated 2 x 10^ times with different realizations of the 
measurement errorsFj The density distributio n of the 
2 X 10^ best-fitting solutions is shown in Figure 10 As 



discussed in Section 4.2, we found that even though the 
mock data had no RM effect at all, there is a clear ten- 
dency to "find" solutions near A = 0° or 180°. This 
should raise a concern about the claimed detection of 
the RM effect with A sa 0° or sal80° with a low signal- 
to-noise ratio. The result of our fitting code applied to 
the actual HARPS data (open circle and dashed con- 

We did not use the MCMC algorithm as it would take to 
long to make chains for 10^ data sets, and because we are only 
interested in the best fitting values of v sin and A for each mock 
data set and not the individual confidence intervals. 




-90 90 
X [deg] 

Fig. 10. — Results for simulated data sets with no RM 
effect. Similar to Figure 9, but this time based on the analysis of 
2 X 10^ simulated data sets with no RM effect but with the same 
time sampling and roughly the same RV precision as the HARPS 
data. The gray shades show the density of the best-fitting values of 
vs'miir and A. The contours enclose 68.3%, 95 %, and 99.73% of the 
best-fitting values. The solid circle shows the [Triaud et alT] | |2010[ l 
result. The open circle and the dashed contours show our results 
of fitting the actual HARPS dataset with our MCMC routine. 



tours in Figure 10) gives values for wsini^, and A that 
are within the area containing 95% of the mock-data so- 
lutions. In this sense the "false alarm" probability (the 
odds of finding such an apparently significant retrograde 
orbit when fitting only random noise) is at least 5%. It is 
probably higher, when one considers that the true noise 
may not be uncorrelated and Gaussian. We therefore 
conclude that the current data do not provide secure in- 
formation on the orientation of the stellar spin relative 
to the orbital spin. 

5. DISCUSSION AND SUMMARY 

We have presented two nondetections of the RM ef- 
fect for the transiting planets WASP-lb and WASP-2b. 
In both cases we gathered high-resolution, high-signal- 
to-noise ratio spectra on nights spanning transits, using 
multiple large telescopes. For WASP-1 there is a weak 
indication of a prograde RM effect, and for WASP-2 we 
did not detect the RM effect. Due to the differences in 
the transit geometry, and in the stellar type, we arrived 
at different conclusions about the relative orientation of 
the stellar spin and orbit in each case. 

Because the transit of WASP-lb has a very low impact 
parameter, the only way to produce a low-amplitude RM 
effect is to have nearly perpendicular sky projections of 
the spin and orbital axes (implying a large misalignment 
in the sky plane), or to have a very low wsinz^,. The 
latter option also implies a likely misalignment, because 
the resulting upper limit on v sin ii, is lower than the ex- 
pected V for a star of the given age and mass. A similar 
comparison can be made between the expected v and the 
lower V sin i^, that is estimated from the breadth of spec- 
tral absorption lines. Thus the data give strong evidence 
for misalignment, although it is not certain whether the 
misalignment is mainly along the line of sight, or in the 
sky plane, or both. 

For WASP-2b, no information on A was gained from 
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our nondetection, mainly because this star is expected 
to be a slow rotator. The upper limit on v sin i^, from 
the RM nondetection is within the expected range of v 
for a star of the given mass and age. An analysis of 
previous HARPS data favored a retrograde orbit for the 
planet, but wc have argued that this may have been a 
statistical false alarm. Numerical experiments confirm 
that fitting random noise with an RM model can produce 
false detections with nearly the same amplitude as the 
claimed detection. For a firmer conclusion one would 
need to gather more spectroscopic data during transits. 
These same numerical experiments should lead to a re- 
evaluation of other cases in which the RM effect was 
det ected with low sta tistical significance, such as TrES- 
2b dWinn et al.|[2008| . 

We now pu t these results into the context of the pat - 
tern noted by Win n et aL] (|2010) and"Schl aufman| ( |2010[ ), 
that hot stars tend to have nigh obliquity, i'he proposed 
boundary line between "hot" and "cool" star was around 

Teff = 6250 K. 

For WASP-2, Cameron et al. (2007) measured an effec- 
tive temperature of 520U ± 2U0 KTand from our HIRES 
spectrum we found 5206 ± 50 K. Thus there is consensus 
that WASP-2 is a cool star. The finding of a retrograde 



eTaL] ( |2007l ), who found T^s = 6110 ± 45 K. Our anal- 
ysis of a HIRES spectrum gave T^g = 6213 ± 51 K, or 
100 K hotter than the determination by Stempels et al. 
j^OTl). Probably th e reason for the difference is that 
[Stempels e t al.|(|2007 ) used the Ha line profile as the main 
constraint on Jeff, while our analysis used th e standard 
SME wavelen gth intervals which exclude Ha (Valenti & 



orbit by |Triaud et al.| ( |2010[ ) was a strong exception to 
the proposed pattern! Uur data and our analysis led 
us to conclude that the spin-orbit angle for this system 
is undetermined, and therefore that WASP-2 is not an 
exception. 



For WASP-1, Cameron et al. (2007) measured an effec- 
tive temperature of 6200 ± 20(J K. Further ob servations 
and spectroscopic analysis were presented by [Stempels 



Fischer 2005 Table 3). It is beyond the scope of this 
article to evaluate the relative merits of these different 
methods for establishing an accurate effective tempera- 
ture scale. Instead we note that the SME-based scale 
that we have used is similar or identical to the scale that 
has been used for the other transit-hosting stars, and 
therefore the scale on which the proposed boundary of 
6250 K is relevant. In this light it seems that WASP-1, 
with Tcff (SME) = 6213 ± 51 K is very near the bound- 
ary. Therefore the finding of a high obliquity neither 
corroborates nor weakens the proposed pattern, although 
WASP-1 may serve as a useful point in establishing the 
sharpness of the transition from mainly-misaligned to 
mainly-aligned. 

We thank G. Marcy and M. Holman for help gath- 
ering some of the data presented here. We are grate- 
ful to the anonymous referee for a prompt and insight- 
ful report, and to Amaury Triaud for comments on the 
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fellowship from the Netherlands Organization for Sci- 
entific Research (NWO). J.N.W. acknowledges support 
from a NASA Origins grant (NNX09AD36G). This re- 
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APPENDIX 



RM effect for WASP-1 and WASP-2 



TABLE 3 

Relative Radial Velocity measurements op WASP-1 



limc [LSJJJtdbJ 


ItV [m s J 


1 

Unc. [m s ] 


Spectrograph 


2454345.83725 


33.86 


2.82 


HIRES 


2454345.84471 


26.32 


2.57 


HIRES 


2454345.85916 


29.23 


2.40 


HIRES 


2454345.86449 


19.53 


2.94 


HIRES 


2454345.86988 


16.26 


2.76 


HIRES 


2454345.87525 


28.59 


2.71 


HIRES 


2454345.88059 


26.28 


2.75 


HIRES 


2454345.88603 


18.84 


2.76 


HIRES 


2454345.89146 


27.23 


2.63 


HIRES 


2454345.89684 


18.61 


3.05 


HIRES 


2454345.90219 


23.13 


2.87 


HIRES 


2454345.90758 


13.82 


2.92 


HIRES 


2454345.91296 


18.70 


2.84 


HIRES 


2454345.91830 


19.45 


2.65 


HIRES 


2454345.92368 


-0.07 


2.78 


HIRES 


2454345.92907 


7.22 


2.77 


HIRES 


2454345.93453 


-8.89 


2.73 


HIRES 


2454345.93990 


2.32 


2.74 


HIRES 


2454345.94841 


1.12 


2.81 


HIRES 


2454345.95376 


-1.71 


2.91 


HIRES 


2454345.95912 


-11.69 


2.60 


HIRES 


2454345.96926 


-5.30 


2.73 


HIRES 


2454345.97461 


-11.05 


2.85 


HIRES 


2454345.97998 


-12.33 


3.32 


HIRES 


2454345.98537 


-10.78 


3.19 


HIRES 


2454345.99861 


-15.13 


2.93 


HIRES 


2454346.00397 


-18.29 


2.98 


HIRES 


2454346.00936 


-12.28 


3.03 


HIRES 


2454346.01468 


-19.91 


2.99 


HIRES 


2454346.02006 


-16.84 


2.99 


HIRES 


2454346.03551 


-22.71 


3.01 


HIRES 


2454346.06736 


-33.51 


2.97 


HIRES 


2454346.13628 


-62.64 


2.96 


HIRES 


2454346.14171 


-64.20 


2.78 


HIRES 


2454318.09458 


63.61 


9.84 


HDS 


2454318.12329 


48.93 


9.86 


HDS 


2454318.13785 


36.74 


10.64 


HDS 


2454350.88327 


51.29 


5.81 


HDS 


2454350.89783 


45.61 


5.81 


HDS 


2454350.90899 


41.88 


6.36 


HDS 


2454350.91661 


31.89 


6.08 


HDS 


2454350.92423 


44.93 


6.19 


HDS 


2454350.93185 


32.72 


6.31 


HDS 


2454350.93945 


40.32 


6.57 


HDS 


2454350.94707 


35.60 


7.06 


HDS 


2454350.95469 


30.06 


5.77 


HDS 








rlUo 


2454350.96992 


32.04 


6.13 


HDS 


2454350.97753 


21.08 


5.92 


HDS 


2454351.03561 


6.90 


9.11 


HDS 


2454351.05624 


-5.33 


7.87 


HDS 


2454351.06385 


-5.00 


7.65 


HDS 


2454351.07146 


-17.63 


8.04 


HDS 


2454351.08247 


-15.36 


6.90 


HDS 


2454351.09703 


-21.98 


6.42 


HDS 


2454351.11159 


-20.25 


8.19 


HDS 


2454351.12615 


-29.56 


7.36 


HDS 



Albrecht et al. 



TABLE 4 

Relative Radial Velocity measurements op WASP-2 



iimc [BJDtdbJ 


T>\ T r^, ... — li 
HV [m s J 


1 

Unc. [m s ] 


Spectrograph 


2454348.72936 


47.57 


6.15 


HDS 


2454348.73875 


46.93 


5.29 


HDS 


2454348.74635 


41.88 


5.75 


HDS 


2454348.75397 


37.82 


5.51 


HDS 


2454348.76158 


23.31 


5.67 


HDS 


2454348.76920 


28.01 


5.53 


HDS 


2454348.77680 


30.28 


5.66 


HDS 


2454348.78442 


17.81 


5.26 


HDS 


2454348.79204 


29.38 


5.52 


HDS 


2454348.79965 


10.77 


5.73 


HDS 


2454348.80728 


6.57 


6.10 


HDS 


2454348.81489 


11.95 


5.42 


HDS 


2454348.82250 


2.32 


5.23 


HDS 


2454348.83352 


-6.93 


5.20 


HDS 


2454348.84809 


-4.37 


4.45 


HDS 


2454348.86264 


-11.82 


4.52 


HDS 


2454348.87720 


-16.39 


4.75 


HDS 


2454348.89175 


-29.35 


4.59 


HDS 


2454348.90632 


-33.89 


4.52 


HDS 


2454348.92089 


-44.93 


4.18 


HDS 


2454348.95000 


-54.24 


4.86 


HDS 


2454363.74812 


79.94 


6.36 


HDS 


2454363.76268 


73.84 


6.94 


HDS 


2454363.79521 


53.93 


7.07 


HDS 


2454363.80283 


50.60 


7.88 


HDS 


2454363.81044 


45.11 


7.11 


HDS 


2454363.81804 


40.19 


6.31 


HDS 


2454363.82566 


41.40 


7.46 


HDS 


2454363.83326 


22.31 


8.60 


HDS 


2454363.84088 


34.76 


7.03 


HDS 


2454363.84850 


28.17 


6.82 


HDS 


2455435.53391 


56.65 


4.38 


PFS 


2455435.54192 


54.05 


4.56 


PFS 


2455435.54978 


45.81 


4.58 


PFS 


2455435.59753 


26.70 


4.62 


PFS 


2455435.60576 


28.47 


6.46 


PFS 


2455435.61657 


28.78 


7.18 


PFS 


2455435.62090 


14.48 


6.44 


PFS 


2455435.62533 


33.21 


6.00 


PFS 


2455435.62971 


13.14 


6.00 


PFS 


2455435.63414 


24.86 


6.09 


PFS 


2455435.63854 


-3.11 


5.62 


PFS 


2455435.64295 


22.52 


5.21 


PFS 


2455435.64731 


12.40 


5.33 


PFS 


2455435.65172 


18.86 


5.71 


PFS 


2455435.65618 


7.36 


5.93 


PFS 


2455435.66062 


-15.67 


5.78 


PFS 


2455435.66496 


0.39 


5.01 


PFS 


2455435.66941 


-12.93 


5.32 


PFS 


2455435.67376 


0.00 


5.24 


PFS 


2455435.67817 


0.19 


5.95 


PFS 


2455435.68261 


1.08 


4.92 


PFS 


2455435.68702 


-27.79 


4.90 


PFS 


2455435.69140 


-4.66 


4.44 


PFS 


2455435.69580 


-6.84 


4.64 


PFS 




1 'X "XGl 
— io.ay 


o.oy 


r^r o 


2455435.70460 


-13.75 


5.15 


PFS 


2455435.70901 


-27.76 


4.88 


PFS 


2455435.71338 


-63.64 


5.72 


PFS 


2455435.71787 


-17.76 


6.67 


PFS 


2455435.72222 


-16.32 


6.81 


PFS 


2455435.72666 


-41.92 


5.87 


PFS 


2455435.73313 


-31.84 


4.42 


PFS 


2455435.74091 


-48.37 


4.57 


PFS 


2455435.74864 


-49.44 


4.97 


PFS 


2455435.75692 


-58.54 


5.51 


PFS 



